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We theoretically consider g factor and spin lifetimes of holes in a longitudinal Ge/Si core/shell 
nanowire quantum dot that is exposed to external magnetic and electric fields. For the ground states, 
we find a large anisotropy of the g factor which is highly tunable by applying electric fields. This 
tunability depends strongly on the direction of the electric field with respect to the magnetic field. 
We calculate the single-phonon hole spin relaxation times Ti for zero and small electric fields and 
propose an optimal setup in which very large Ti of the order of tens of milliseconds can be reached. 
Increasing the relative shell thickness or the longitudinal confinement length further prolongs T\. 
In the absence of electric fields, the dephasing vanishes and the decoherence time T2 is determined 
by T 2 = 2Ti. 

PACS numbers: 71.70.Ej, 63.22.-m, 81.07.Vb, 81.07.Ta 



Semiconducting nanowires (NWs) allow to create 
nanoscale systems denned precisely regarding composi- 
tion, geometry, and electronic properties and hence are 
subject to great experimental efforts. Furthermore, they 
offer new ways for implementing spin-based quantum 
computation. 1 Both III-V compounds and group-IV ma- 
terials are considered and operated in the conduction 
band (CB, electrons and in the valence band (VB, 
holes)ME2 re gi me A particularly favored material is 
In As, where single-electron quantu m do ts (QDs)^ and 
electrically controlled spin rotation d 5 * 6 * 8 ! have been im- 
plemented. Recently, qubits have also been implemented 
in InSb NW QDs,-22I2! a system for which extremely 
large electron g factors have been foundP^ However, the 
strong hyperfine interaction in InAs and InSb is consid- 
ered the dominant source for the short coherence times 
observedP^ The latter may therefore be substantially 
prolonged in group-IV NWs that can be grown nuclear- 
spin-free. In this context, Ge and Si have emerged as 
promising materials for nanoscal e sy stems such as lateral 
QDsjMM] S elf-assembled QDsj22H22l cylindrical core/shell 
j^YgCflHio] and ultratmri) triangular NWsP 

For applications in spintronics and quantum informa- 
tion processing, it can be advantageous to consider holes 
instead of electrons. Due to the p-wave symmetry of the 
Bloch states, holes experience a strong spin-orbit inter- 
action (SOI) on the atomic level leading to an effective 
spin J = 3/2 behavior. Hence spin and momentum are 
coupled strongly which allows efficient control of the hole 
spin by electrical means. Furthermore, hole spin lifetimes 
are prolonged in the presence of confinement OSSHU] 

In Ge/Si core/shell NWs, the la rge V B offset leads to 
an accumulation of holes in the core EEU They form a one 
dimensional (ID) hole gas with an unusually large, tun- 
able Rashba-type SOI, referred to as direct Rashba SOI 
(DRSOI)P This DRSOI makes Ge/Si core/shell NWs 
attractive candidates for quantum information process- 
ing via electric-dipole induced spin resonance,^ and we 
mention that signatures of a tunable Rashba SOI were 
already deduced from magnetotransport experiments.^ 



Experiments on gate defined QDs in this system re- 
vealed a n anis otropy and confinement dependence of the 
g factorP^IISl Recently, singlet-triplet relaxation times 
in the range of several hundred microseconds were 
measured.^ 

In this work, we consider holes forming qubits in the 
energetically lowest states of longitudinal QDs in Ge/Si 
core/shell NWs. We find the effective g factor g e g of this 
subsystem which turns out to be strongly anisotropic and 
tunable by choosing the direction and magnitude of ap- 
plied electric fields. For small electric fields, we perturba- 
tively derive an effective subspace Hamiltonian and the 
according hole spin phonon coupling and calculate the 
hole spin relaxation rate T-j - . At small Zeeman split- 
tings hu we observe a w 7 / 2 proportionality of T^ 1 which 
contrasts the oj 5 behavior found for electrons in QDsPSES] 
The magnitude of T^ 1 depends strongly on the direction 
of the magnetic field with respect to the wire. For zero 
electric field, aligning the magnetic field perpendicular 
to the wire results in very long Ti of the order of tens 
of milliseconds. Directing the magnetic field along the 
wire results in a much shorter T\. For both configura- 
tions, the dephasing is zero, hence the decoherence time 
is given by T 2 = 2T\. Applying small electric fields can 
enhance the relaxation rate by several orders of magni- 
tude. This effect depends strongly on the direction of the 
electric field with respect to the magnetic field. Long Ti 
in the presence of electric fields are obtained when elec- 
tric and magnetic fields are perpendicular to each other 
and perpendicular to the wire. Moreover, we find that Ti 
can be prolonged further by increasing the relative shell 
thickness and the longitudinal QD confinement. Thus, 
we predict an optimal field geometry for spin qubits in 
Ge/Si NWs that can be tested experimentally. 

Low-energetic hole states in a cylindrical Ge/Si 
core/shell NW are well described by an effective ID 
HamiltoniarP^ 

H w =H a + H' (1) 
that can be split into a leading order term Hq and a 



perturbation H', 



Ho — i?LK d 
H = J?LK„ 



H strain + Hb,Z, (2) 
I" + i/DR + #3, orb- (3) 



Using the notation introduced in Ref. 3? and denning the 
z axis as the NW axis (see Fig. [I]) , the diagonal terms 
of the Luttinger-Kohn (LK) Hamiltonian and the strain- 
induced energy splitting read 



HhK d + -Hsi 



A+(k z ,j) + A-(k z ,-y)T z 



(4) 



Here, tj and Oi are the Pauli matrices for band in- 
dex ({g,e}) and spin block ({+,—}) of the basis states 
9±( x jV) an d e±(x,y) that provide the transverse mo- 
tion. In Eq. pi, we defined A±(k z ,j) = f?k 2 z {mZ x ± 



1 ±l'"zi I) — '"zV""j 

m ( T 1 )/4 ± A/2, with m g ~ Wo/(7i + 27 s ) and m e = 
"^0/(71 +7s) as the effective masses along z. Here, 71 and 
7 S are the Luttinger parameters in spherical approxima- 
tion and mo denotes the bare electron mass. For Ge, 71 = 
13.35 and 7 S = 5. II. 45 A = Alk + A st rain(7) is the level 
splitting between the g± and e± states, 7 = (R s — R)/R is 
the relative shell thickness, and R (R s ) is the core (shell) 
radius. The Zeeman coupling Hb,z with splitting hwB.z 
in the lowest-energy subspace (g band) is determined by 
the magnetic field B = (B x ,0,B z ) = \B\(sin9, 0, cos 6) 
(Fig. [T]), where we set B y = due to cylindrical symme- 
y. the 



try. the main contributions to H' are 



HhK od = Ck z T y a x , 
H D r = eU(E x T x a z 



EyTy), 



(5) 
(6) 



where i?LK od features the off-diagonal couplings with 
coupling constant C — 7.26h 2 / (m R) provided by the LK 
Hamiltonian as a consequence of the strong atomic level 
SOI. .Hdr is the DRSOI that results from direct, dipo- 
lar coupling to an electric field E = (E x , E y , 0), where 
U = 0.15R. We note that hk z = -ihd z in Eqs. Q and 
([5]) is the momentum operator along the wire. In the 
absence of longitudinal confinement the wave functions 
along z are of type e lk * z with k z as the wavenumber. 
H-r is the conventional Rashba SOI, and, although fully 
taken into account in the present analysis, turns out to 
be negligible for the typical parameters and electric fields 
considered here. Finally, Hb, orb denotes the orbital cou- 
pling to the magnetic field. Details on all ele men t s of Hp 
and H' are provided in Ref. 37 and in Eqs. (A1|-(A7| in 
the appendix. 

We proceed with the derivation of an effective 
ID Hamiltonian i/h-ph for the coupling between low- 
energetic holes and acoustic phonons. There are three 
different types of acoustic phonon modes in cylindric 
NWs: torsional, dilatational, and flexuralP^ We find four 
different modes A with dispersion relation oj\(q), where 
q is the phonon wavenumber along the wire and the 
exact form of u\ depends strongly on the shell thick- 
ness. For the torsional and dilatational mode (A = T, L) 
oj\ depends linearly on q, whereas for the two flexural 
modes (A = F±i) this dependence is quadratic. The 
detailed derivation will be published elsewhere, in this 




Figure 1. Sketch of a Ge/Si core/shell NW aligned with the 
z axis of the coordinate system. Electric gates (blue) induce 
confinement along the z axis and define a QD. The electric 
field E lies perpendicular to the wire in the xy plane and the 
magnetic field B lies in the xz plane. 



work we directly apply the displacement field u(r,t) = 
J2\ q [ u \{<li r : t)bq,\(t) + H.c] obtained for a finite shell 
following Refs.lMgi Here, b q , x (t) = e~ lul ^ q H q ^ is the 
time-dependent phonon annihilation operator. To derive 
^h-ph, we insert the associated strain tensor components 
£ij(r,t) in the Bir-Pikus Hamiltonianp^l 



H B p=b 



^2 £ H J i + 2 ( £ xy{ Jx, Jy} + Cp.) 



(7) 



where we omitted the global shift in energy and used 
the spherical approximation. The J i? i = x,y,z, are the 
effective spin-3/2 operators of the VB electrons and the 
anti-commutator is defined as {A,B} = (AB + BA)/2. 
For Ge, the deformation potential b takes the value b ~ 
-2.5 eVPWe finally obtain 



h-ph 



5> 



Hi 



Hi 



(8) 



by integrating out the transverse part of the matrix ele- 
ments, i.e. by projecting the Hamiltonian onto the sub- 
space spanned by g± and e± . The components of i?h_ p h 
are given explicitly in Eqs. (B1)-(B4| in the appendix. 



Longitudinal confinement is realized by electric gating 
(see Fig. [TJ , which is modeled by adding a harmonic con- 
finement potential in the z direction, 



H, 



qd 



Vc(z), 



(9) 



where V c (z) — ^a c z 2 . H q d describes the QD well if the 
longitudinal confinement length is much larger than R. 
The basis states of H qc [ are products of type g±il'm an d 
e ±V'm' wnere the ipm £ (z) are eigenf unctions of the har- 
monic oscillator h 2 k 2 z / {2m g / e ) + V c (z) and m e {0, 1, . . .} 
is the harmonic oscillator quantum number. The con- 



finement energies hw g / e relate to a c via a c — 



l 9/e UJ g/ e 



and the harmonic oscillator confinement lengths read 

z g/e = \/fl/(m g / e UJ g/e ). 



3 




Figure 2. Effective g factor g e g as a function of the angle 
defined by B = |B|(sin0, 0, cos 0) for (a) and E\\y 

(b). We vary \E\ from to fO V/fim. For |B| = 0, we 
find <7cff(0) ~ 0.f4 and g c s(-K/2) « 5.7. ft is clearly visible 
that <7 c ff is affected much stronger by changes in \E\ for 
than for E\\y. Even though the curves in (a) seem to overlap 
for \E\ > 6 V//im, g c g(ir/2) still decreases for growing fields 
and g ff remains anisotropic. We choose R = 10 nm and 
R s = f 3 nm for the NW and a QD confinement length of 
z„ « 80 nm. 



From H q d we extract the effective g factor g e ff of the 
lowest-energy subsystem by performing an exact, nu- 
merical diagonalization which gives the Zeeman splitting 
A-Ez,num (defined as positive) and 



AE 



Z,num 



Hb\B\ 



(10) 



where [Ib denotes the Bohr magneton. In Fig. [2j we plot 
g e g as a function of the angle 9, for both E\\x and E\\y. 
In both cases, g e g is highly anisotropic and tunable over 
a wide range of values by adjusting the magnitude of E. 
The tunability is caused by two mechanisms which occur 
in the system for large \E\. The admixture of the e± 
states to the effective lowest-energy subsystem increases 
while the spin-orbit length Zgoi decreases. For very small 
'soi {hoi "C Zg), the hole spin flips many times while 
moving through the QD and the resulting g e g starts to 
average out. The tunability is much stronger for E\\x 
than for E\\y. Note that g e g is also tunable by varying 
V c (z). We find good agreement with the results given in 
Ref. HH where ffoxp ~ 1-02 was measured for B aligned 
with the NW with an accuracy of ~ 30°. We note, how- 
ever, that clearly different results for g can be expected 
in QDs with very large occupation number, i.e. when the 
hole spin qubits are formed in an excited band. 

In the following, we are interested in the dynamics of 
the lowest lying, Zeeman split states which we decou- 
ple perturbatively from the higher energy states. This 
is done by two consecutive Schricffcr- Wolff transforma- 
tions (SWTs) to account for the two different energy 
scales A and hui g . The general form of the SWT is 
H = e~ s He s , where to lowest order S ~ Si. We first 
remove the coupling between the g± and e± states in 
the effective ID picture using Sf. The hole-phonon cou- 



pling then transforms according to i/h-ph — [Sf , -ffh-ph] 
and we refer to its projection on g± as H^_ ph . In the 
second step, we add harmonic confinement as introduced 
above and decouple the two lowest, Zeeman split states 
|0) = {|fr), HJ.)} by another SWT using s[ 0) . A necessary 
condition for this approach is that the energy splittings 
obey A 3> hu> g ^> Hujb,z, and the magnitude of E is re- 
stricted by 2C\E\eU / (z g A) <C huj g . The latter condition 
is fulfilled for \E\ <C 1 V///m. We obtain an effective 
Zeeman term Hz, c g = /is-B c ff ■ cr with Zeeman splitting 
AEz,eS = 2/i_B|-B c ff|; where a is a vector of Pauli matri- 
ces. The effective hole spin phonon coupling is obtained 
by taking 



ph 



H 9 

-"h-ph 



[S[ 



10} 



H 9 1 



(11) 



where ph is now written in the basis given by the 
confinement. Projecting -ff s _ p h on |0) results in an effec- 
tive coupling -ff s -ph,off = /J-b$B ■ cr with the fluctuating 
magnetic field 8B(t) = J2x. q [a\(q)b q ,\(t) + H,.c.]. The 
effective subspace Hamiltonian then reads 



H ff = H 



X.eff + -ffs- P h,cff = Ms {B c g + 6B(t)) • cr. (12) 

The spin relaxation rate in the Born-Markov approxima- 
tion is given by the Bloch-Redfield approac h 1 40 * 50 * 51 ! 

1 

Ti~ : 



[Sij(S pq - n p n q )J+(uj) - (6 ip 



-Sij£kpqn k Ip g (u) + E ipq n p I qj (ui)] , 



(13) 



where summation over repeated indices is assumed, n = 
B e g/\B e g\ is the unit vector in direction of the effective 
magnetic field, and hu = hioz. c g = AEz, e fi is the en- 
ergy splitting of the considered states. Here, J^(uj) = 
Re[Jy(w) + Jij(-uj)] and I^{uj) = Im[Jy(a;) - Jy(-o;)], 

with Jij(u;) = {^b/K) 2 r o °° dte-^iSBi^SBjit)) denot- 
ing the spectral function. 

In Fig. [3J we display T^ 1 for \E\ —Q and two different 
directions of B with respect to the wire, B\\z and B\\x. 
In this case, the spin-phonon coupling -ff s -ph,eff depends 
only on the coupling terms of i?s j0 rb- For low u>z, c s, 
i.e. the long wavelength regime (qz g <C 1), both curves 
7/2 

are proportional to uj z eS . This behavior is valid for low 
temperatures (Hujz,cS ^S> ksT) and will be replaced by 
Tj~ oc u> z eS T for hu>z,cS "C ksT. The ui ZeS scaling con- 
trasts the lo z eff behavior of electrons in QDsPSES p or 

B\\z, only the F±i modes contribute significantly to . 
When directing B||£C, the F±i contributions dominate 
for low luz.cS and, for the chosen QD geometry, are re- 
placed by a dominating T contribution at \B\ w 150 mT 
(Fig. [3j dashed) . This results in a double peak whose rel- 
ative height can be modified by changing z g or R and R s . 
Most remarkably, for B\\x, T^ 1 is several orders of mag- 
nitude smaller than for B\\z. For the chosen QD geom- 
etry, T-f 1 reaches maximal values T^^J^BWx) ps 60 s _1 

and rj~ max (B||i) rj 11 ms _1 . These rates are, depend- 
ing on the direction of B, comparable to or much smaller 
than for electrons in In As NW QDs.^ 
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Figure 3. Relaxation rate T^ 1 at \E\ = for B\\x (red, 
solid) and B\\z (blue, solid). For -B||a;, we plot the contribut- 
ing phonon branches F±i and T (dashed). We find maximal 
values T-^(B\\z) « 11 ms" 1 and r^(B||x) » 60 s" 1 . 
Note the non-monotonic behavior of T-j -1 as function of uiz,eB- 
The NW and QD parameters are chosen as in Fig. [2] 
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Figure 4. Relaxation rates T{~ for E\\x with \E\ — 0.1 V//xm 
for B||& (red, solid) and B\\z (blue, solid). For comparison 
we plot Tj -1 at \E\ — (dotted). We find maximal values 
^LfBp) « 3.2 ns~\ ^(Bp) a 5.8 ms" 1 . Rotating 
the electric field so that E\\y yields the same curve for B\\z. 
Remarkably, for B\\x almost no difference between the curves 
at E\\y (dashed) and \E\ = (dotted) is observed. We use 
the NW and QD parameters given below Fig. [2] 



larger u>z,cS and lowers the peak height, e.g. increasing 7 
from 0.3 to 0.7 reduces by a factor ~ 3. However, 
decreasing (increasing) R and R s while keeping 7 con- 
stant has no substantial effect on aside from slight 
shifts to the right (left) on the u>z, e s axis. Additionally, 
for \E\ = and for E\\y and B\\x, enhancing the confine- 
ment length z g lowers T^^^ since the short wavelength 
regime is reached for smaller (Jz,eB- This effect is quite 
large, for instance raising z g from 60 to 100 nm tunes 
^lTmax by factors between 10 and 100. However, for E\\x 
and B\\z or B\\x, increasing z g results in larger T^~^ ayi . 
From this analysis we conclude that there exist optimal 
configurations of B and E in order to obtain long T\ in 
this type of NW QD. B should be applied perpendicular 
to the NW and the optional E should lie perpendicular to 
both B and the NW. For vanishing B, as pointed out in 
Ref. two-phonon processes^ might become relevant. 

In the Bloch-Redfield framework, the decoherence time 
is given by T" 1 = (2 T1)- 1 + T~ x , where T v denotes the 
dephasing time.^MJ For \E\ = 0, we find T 2 = 2T X be- 
cause the corresponding spectral function is super-ohmic 
and gives T" 1 = 0. For \E\ ^ the SOI results in a 
non-zero dephasing term T" 1 7^ and hence T 2 < 2T X . 

In conclusion, we have examined effective Zeeman 
splitting and hole spin dynamics for holes in the lowest 
VB of a Ge/Si core/shell NW QD. We reported a highly 
anisotropic effective g factor which is strongly tunable 
by applying electric fields. We calculated relaxation 
rates and found configurations of electric and magnetic 
fields which correspond to very long spin relaxation 
times. Furthermore we pointed out that the relative 
shell thickness and the QD confinement length influence 
the spin relaxation time. 

We thank Peter Stano for helpful discussions. This 
work has been supported by SNF, NCCR Nano, NCCR 
QSIT, DARPA, and IARPA. 

Appendix A: Effective ID Hamiltonian 



Considering non-zero electric fields, we plot T-f 1 for 



We add the cor- 
dashed) to allow for 
B, T^ 1 is enhanced 



E\\x again for B\\z and B\\x (Fig. 
responding curves for \E\ =0 (Fig. [4 
comparison. For both orientations o:' . 
significantly for larger wz.eff- This is due to phonons of 
the L mode coupling |fr) and |JJ.) via a combination of 
i?LK od and i?DR. which dominates Hr. Due to cylin- 
drical symmetry, applying E\\y for B\\z results in the 
same effect as described for E\\x. Remarkably, in stark 
contrast to £?||J5||i, only minor changes with respect to 
the curve at \E\ — (Fig. |4j dotted) are observed when 
E\\y and B\\x (Fig. Hi dashed). In the latter case, the 
dominant contributions of -ff s _ p h.cff are already present in 
#s- P h,eff for \E\=0. 

In both cases, \E\ — and \E\ 7^ 0, increasing the 
relative shell thickness 7 shifts the curves to slightly 



Here we display the effective ID Hamiltonians used in 
Eqs. ^ and ([3f of the main text. An extended deriva- 
tion of these terms can be found in Ref. [371 We use the 
basis {3+, 3-, e+, e_}, where the exact form of the states 
9±i x iy) and e±(x,y) is given in Ref. l37l The diagonal 
elements of the Luttinger-Kohn (LK) Hamiltonian and 
the strain induced splitting is combined in 



H 



LK d 



strain — 





h 2 k 2 z 




2m g 












V 








h 2 kl 

2m g 






hH-j 

2m, 



2m t 



A 



. (A1) 

where hk z is the momentum operator along the wire and 
frig ~ mo/ (71 + 27 s ) and m e = mo/ (71 + J s ) are the ef- 
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fective masses along z. Here, 71 and 7 S are the Luttinger 
parameters in spherical approximation and mo denotes 
the bare electron mass. A = A LK + A strain (7) denotes 
the level splitting between g± and e±, 7 = (R s — R)/R is 
the relative shell thickness, and R (R s ) is the core (shell) 
radius. Alk = 0-73H 2 /(moR 2 ) and the strain dependent 
energy splitting can take values A stra ; n (7) ~ — 30 meV. 
The magnetic field B = (B x ,0,B z ) = |B|(sin0, 0, cos 0) 
interacts via the Zeeman coupling 



Hb,z — Ms 



/ B Z G 
B X K 


V 



B X K 
~B Z G 










B Z F B X M 
B T M —B,F 



(A2) 



where we set B y = due to cylindrical symmetry. Here, 
/is is the Bohr magneton and the parameters take the 
numerical values G = -0.06, K = 2.89, M = 2.56 and 
F = 1.56. From the LK Hamiltonian additional off- 
diagonal coupling terms arise and we find 



LK od 







\iCk z 






iCk z 






-iCk z 






(A3) 



pling to E is given by 



Hur = eU 



( 





E 


\ 











-E* 


E* 











V 


—E 





/ 



(A5) 



where the parameter U = 0.15i? scales linearly in R. We 
note that the parameters S and U of Hr and Hur are 
related by 



eU 



-1.1 



R 2 



(A6) 



hence i?DR dominates by one to two orders of mag- 
nitude for radii between R = 5 — 10 nm. Lastly, the 
interaction of the magnetic field via the orbital motion is 
given by 

/ -B x Lk z ~B z Dk z 

_ B z Dk z -B x Lk z 

UB ' OTh ~ B x Lk z -B z Dk z 

\ B z Dk z B x Lk z 



(A7) 



with L 



MR and D = 2.38R. 



Appendix B: Effective ID phonon Hamiltonian 



with coupling constant C = 7.26H 2 /(m R). To describe 
interactions with the electric field E — (E x ,E y ,0) — 
(cos t/7 c i, sin i^ei, 0) , we introduce the shorthand nota- 
tion E = |£J|e 4< ^ cl . The effective conventional Rashba 
SOI interaction reads 



a 



/ 

iE*k z T 
E*S 

V 



-iEk z T 



-ES 



ES 



\%Ek z 





-E*S 
■\iE*k z 




(A4) 



with T = 0.98 and S = 0.36/R and a = -0.4 nm 2 e. 
Although fully taken into account, turns out to be 
negligible for the typical parameters and electric fields 
considered here (see Eq. (A6|). The direct, dipolar cou- 



Starting from the spherical Bir-Pikus Hamiltonian (see 
Eq. (7]) of the main text), we derive an effective hole 
phonon coupling Hamiltonian for each phonon mode A. 
This is done by integrating out the transverse compo- 
nents of the matrix elements of the states g±{x,y) and 
e ±{ x ilj)- To improve readability, we introduce an ef- 
fective phonon annihilation operator b\ — e lqz b q x(t) — 

The transversal phonon mode T couples as 



H 7 



E 



01 






b T -U 



-bn 



and the dilatational mode L gives 







(Bl) 



Hr 



E 



( a 2 (b L + b{) 

a 2 (b L + b\) ia 3 (b L 

ia 3 (i 



V ia 3 (b L 



St 



% 

A) 



a 4 (6 L + b 







ia 3 (b L -b{)\ 

~A) 

1) 

a 4 (b L +b[) J 



(B2) 



For the two flexural modes F ±1 , we find 



H,, 



1 



ia 5 






~b F+ 



\ 









\ 



(B3) 



and 



H, 



E ia s 



( 





-6f_! 


























V 


b F _ 1 





/ 



(B4) 



G 



Here, the aj, i = 1, 2, 3, 4, 5, are real, q-dependent prefac- 
tors. The complete effective hole phonon coupling ifh-ph 
is then given by 

-ffh-ph = ^ H\ = Ht + i?L + Hp +1 + Hf_! ■ (B5) 

A 
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